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Abstract 

We show that the probability densities of accelerations of La- 
grangian test particles in turbulent flows as measured by Bodenschatz 
et al. [Nature 409, 1017 (2001)] are in excellent agreement with the 
predictions of a stochastic model introduced in [C. Beck, PRL 87, 
180601 (2001)] if the fluctuating friction parameter is assumed to be 
log-normally distributed. In a generalized statistical mechanics set- 
ting, this corresponds to a superstatistics of log-normal type. We 
analytically evaluate all hyperflatness factors for this model and ob- 
tain a flatness prediction in good agreement with the experimental 
data. There is also good agreement with DNS data of Gotoh et al. 
We relate the model to a generalized Sawford model with fluctuat- 
ing parameters, and discuss a possible universality of the small-scale 
statistics. 



In the past two years there has been considerable experimental progress 
in mesasuring the statistics of Lagrangian test particles in turbulent flows 
|2|, |3|]. While simple dynamical models like the Sawford model ||, [J predict 
a Gaussian acceleration statistics, the measurements of Bodenschatz et al. 
HI, || have confirmed that the probability distribution of acceleration a is 
strongly non-Gaussian and has strongly pronounced tails. As shown in || |7| . 
the measured acceleration distributions are well fitted by Tsallis distributions 
|| for \a\ < 30, in units of the standard deviation. However, very new 
recent measurements |§ and numerical simulations |L0| indicate that there 



are deviations from Tsallis statistics for extremely large accelerations \a\ > 
30, i.e. for extremely rare events of very strong forces acting on the particle. 

In this Letter we show that the simple dynamical model previously intro- 
duced in Jf| can be generalized in such a way that it yields quite a perfect 
description of the observed Lagrangian acceleration statistics in turbulent 
flows. The coincidence is better than in previous models based on Tsallis 
statistics. In the language of statistical mechanics, the relevant new statis- 
tics that fits the data quite perfectly is 'superstatistics of log-normal type'. 
While it can be proved that log-normal superstatistics, as any superstatis- 
tics, is close to Tsallis statistics for not too large accelerations, there are 
significant differences in the tails. 



The concept of superstatistics was introduced in JTT]. A superstatistics 
is a 'statistics of a statistics', one given by ordinary Boltzmann factors e~^ E 
and the other one by the distribution function of (3, which is assumed to 
be fluctuating. Superstatistics are generally relevant for an effective descrip- 
tion of driven nonequilibrium systems with stationary states which possess 
a fluctuating intensive parameter (3. This can, for example, be the inverse 
temperature or, in the turbulence application, the energy dissipation rate. 
Tsallis statistics is a special superstatistics obtained for a ^-distribution 
of the intensive parameter. But as shown by Tsallis and Souza [12], more 
general versions of statistical mechanics can also be constructed for other 
superstatistics than Tsallis statistics, in particular for the superstatistics of 
log-normal type which turns out to be relevant here. 

A superstatistics can be easily dynamically realized, by generalizing the 
approach of [0. Consider a Langevin of the form 

a = jF(a) + cxL(t), (1) 

where L(t) is Gaussian white noise, 7 > is a friction constant, a describes 
the strength of the noise, and F(a) = -£V{a) is a drift force. For most 
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applications, it is sufficient to consider linear forces F(a) = —a as generated 
by V(a) = |a 2 . For an ordinary Brownian particle, a would be the velocity, 
but for the Lagrangian turbulence application, a actually stands for the ac- 
celeration of the Lagrangian test particle, i.e. a velocity difference on a very 
small time scale. We will consider linear Langevin equations of the form ([!]) 
for each component of the acceleration. 

If 7 and a are constant then the stationary probability density of a is 
p(a) ~ e~^ v<yCi \ where (3 := However, since the energy dissipation fluc- 
tuates in a turbulent flow we now let the parameters 7 and a fluctuate such 
that :— -3y has some probability density f((3)- These parameter fluctua- 
tions are assumed to be on a relatively long time scale (or a relatively large 
spatial scale) so that the system can temporarily reach local equilibrium. The 
Lagrangian test particle moves through these different regions with different 
(3. 

One obtains for the conditional probability p(a\0) (i.e. the probability of 
a given some value of (3) 

p(a\(3) = -^e W {-(3V(a)}, (2) 

where Z(/3) is a normalization constant, and for the joint probability p(a, (3) 
(i.e. the probability to observe both a certain value of a and a certain value 
of/3) 

p(a,/3)=p(a\P)f(/3), (3) 

and finally for the marginal probability p(a) (i.e. the probability to observe 
a certain value of a no matter what (3 is) 

p(a) = [ p(a\(3)f((3)d(3. (4) 



In a ^-distribution W as chosen for f((3), and it was shown that this 
generates Tsallis statistics for the marginal distribution, i.e. p(a) ~ (1 + 
$(q— l)V(a))^ 1 ^ q ^ 1 \ where q is the entropic index of nonextensive statistical 
mechanics |§. 

Let us here choose another distribution for f((3), in fact one of the exam- 



ples dealt with in [11], the log-normal distribution 
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where m and s are parameters. The average (3q of the above log-normal 
distribution is given by j3 = rriy/w and the variance by a 2 = m 2 w(w — 1), 

2 

where w := e s . For our Lagrangian turbulence model it is sufficient to 
consider linear drift forces F(a) = —a. The integral given by (|4]) 

p(a) = -L [°° df3 /3- 1 ' 2 exp { ~ (1 " g 2 ^ )2 \ e-& 3 , (6) 



2ns J I 2s 2 

cannot be evaluated in closed form, but the equation is easily numerically 
integrated, and can be compared with experimentally measured densities 



p(a). Similar distributions as in (|6]) have also been considered in | I3| , however 
without a dynamical realization in terms of a stochastic differential equation. 
The distribution p(a) has variance 1 for the choice m = y/w, hence only one 
fitting parameter s remains if one compares with experimental data sets that 
have variance 1. 

Fig. 1 shows the histograms of accelerations of Lagrangian test particles 
as measured by Bodenschatz et al. |], §]. The distributions are rescaled to 
variance 1. As shown previously |], 0], the experimental distributions are 
reasonably well approximated by Tsallis distributions for \a\ < 30. But log- 
normal superstatistics yields a better fit, using only one fitting parameter 
s. The measured distributions for R\ = 690 and R\ = 970 are basically 
the same and very well fitted by eq. (Bf) with s 2 = 3.0. Note that m is not 

a free fitting parameter but fixed as m = y/w = to give variance 1. 
Since Bodenschatz's data reach rather large accelerations a (in units of the 
standard deviation), the measured tails of the distributions allow for proper 
distinction between various superstatistics. The main difference between the 
X 2 -superstatistics studied in and log-normal superstatistics studied here 
is the fact that p(a) decays with a power law for the former ones, whereas 
it decays with a more complicated logarithmic decay for the latter ones. For 



alternative fitting attempts based on multifractal models, see |TJ| . 

All moments of the distribution @ exist and can be easily evaluated as 
follows. The moments of a Gaussian distribution 

_-i/3a 2 




PG(a) = V2^ e (?) 
are given by 

(* r h = ^(r-l)» (8) 
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(r even). Moreover, the moments of the lognormal distribution (|^) are given 
by 

((3 r ) L N = m r w^ 2 . (9) 

Combining eq. @ and (H) one obtains the moments of the distribution p(a) 
given by eq. (0) as 

(a r ) = ((a r ) G ) LN (10) 
= (r-l)\\(p-^) LN (11) 
= (r - l)\\m-^w^ r2 (12) 



In particular, the variance is given by 

(a 2 ) = m~ l ^w. (13) 
All hyperflatness factors F r are independendent of m and given by 

F r := ^1 = (2r - l)!!^" 1 ). (14) 

{a?) r v ' K ' 



In particular, the flatness F 2 is given by 



(« 2 ) 



3iw = 3e s . (15) 



Measuring the flatness F 2 of some experimental data thus provides a very 
simple method to determine the fitting parameter s of lognormal superstatis- 
tics. Or, if s is fitted from the shape of the densities, then a prediction on 
the flatness factor can be given. Note that w is the analogue of q in the 



nonextensive approach [ 1 1 



Using the parameter s = 3.0 that yields a good fit of the densities mea- 
sured by Bodenschatz at R\ = 690 and R\ = 970, a flatness of F 2 = 60.3 
is predicted from eq. (]T5|), in agreement with the measured flatness values 
reported in || [|. 

In direct numerical simulations (DNS) of the Navier-Stokes equations, 
even larger accelerations a can be reached. Fig. 2 shows Gotoh's results 
on the pressure distribution as obtained by DNS at R\ = 380 [|H]]. In good 



approximation the pressure statistics coincides with the acceleration statistics 
of a Lagrangian test particle. Gotoh's histograms reach accelerations up to 
150 (in units of the standard deviation), a much larger statistics than can 
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be presently reached in Bodenschatz's experiment. Hence the tails of these 
distributions can very sensitively distinguish between various superstatistics 
models. Fig. 2 shows that log-normal superstatistics with s 2 = 3.0 again 
yields a good fit, keeping in mind that one compares data that vary over 12 
orders of magnitude. The fit quality can be slightly further improved if one 
uses for f((3) a log-normal distribution that has a lower and upper cutoff 
(Anm ~ 0.0005 and /3 max ~ 25). The above truncation may effectively model 
finite Reynolds number and finite size effects, which are certainly present in 
any numerical simulation of the Navier-Stokes equation. 

If the Reynolds number is sufficiently large (R\ > 300), we notice that 
the same fitting parameter s 2 = 3.0 yields a good fit of all 3 different data 
sets (R\ = 380,690,970). It seems that s 2 is essentially independent of the 
Reynolds number for large R\, in agreement with the measured approximate 
Reynolds number independence of the flatness reported in M for R\ > 300. 
If this tendency is confirmed by further experiments then it seems natural 
to conjecture that the variance of the log-normal distribution converges to 
a finite universal value for R\ — > oo. This value might simply be given by 
s 2 = 3. The value 3 could stand in connection with the three spatial degrees 
of freedom: As it is apparent from eq. (||), the variable Y := log^ is a 
Gaussian random variable with average and variance s 2 . The fluctuating 
(3 is related to the fluctuating energy dissipation rate. Since energy can flow 
independently into the three space directions, we may write 

Y = log^ + log^ + log^, (16) 
m m m 

where each of the independent random variables Xj := log—, % = x,y,z, is 
Gaussian with average 0. Our fitting observation s 2 = 3 is thus equivalent 
to the fact that each of the Gaussian random variables X{, describing the 
energy flow in the i-th direction, has variance 1. Whether this is a random 
coincidence for finite Reynolds number or whether this is a deep fundamental 
principle underlying the small-scale statistics for R\ — > oo is not clear at the 
moment. 

What remains to be done is to precisely relate the parameter /3 of our su- 
perstatistics approach to the fluctuating energy dissipation e in the turbulent 
flow. For this we may just compare with a previously studied Lagrangian 
model, the Sawford model ||, |], [T5[| . The Sawford model assumes that the 
joint stochastic process (a(t), u(t), x(t)) of acceleration, velocity and position 
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of a Lagrangian test particle obeys the stochastic differential equation 



-(T^ + t-^a-T^t-'u 



+^2al(T^ + t^)T£H^L(t) (17) 

u = a (18) 
x = u, (19) 

where L(t) is again Gaussian white noise. T L and t v are two time scales, with 
Tl >> t, q . In this model one has 

2cr 2 

TL = cfe (20) 

tv ~ a^ 2 "' (21) 

(22) 

where e is the average energy dissipation, C ,a are Lagrangian structure 
function constants, and a\ is the variance of the velocity distribution. The 
Taylor scale Reynolds number is given by 



Rx = ^S*. (23) 



The Sawford model predicts Gaussian stationary distributions for a and u, 
and is thus at variance with the recent measurements. However, a straight- 
forward idea is to generalize the Sawford model with constant parameters to 
a generalized Sawford model with fluctuating parameters, following a simi- 
lar type of arguments as in [|7J. This was recently worked out by Reynolds 



||15|| , who emphasized compatibility of the obtained stochastic model with 
Kolmogorov's refined similarity hypothesis. 

We note that there are several possibilities to extend the Sawford model 
with constant parameters to an extended one with fluctuating ones, depend- 
ing on which of the variables e in the above model equations are replaced 
by a fluctuating e and which are not. For our purposes it is sufficient to 
consider the Sawford model in the limit Ti — > oo, which is a reasonable 
approximation for large Reynolds numbers. In that limit the model with 
constant parameters becomes identical to eq. (jj) with F(a) = —a and with 
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the identification 

7 = 9l. v -V*$l* (24) 
2a 

r 3/2 

' = ir"" 1 "* (25 » 

This clarifies the physical meaning of the average of our parameters 7 and a. 
To proceed to a model with fluctuting parameters, the simplest possibility is 
to replace in both of the above equations e by a fluctuating e. This yields 

P=^ = ^ /2 e- 3 ^ (26) 

i.e. the fluctuating energy dissipation e is proportional to /3~ 2 / 3 . Based on 
this equation, one can derive a couple of interesting relations. For example, 
for the constant a$ defined by the relation 

(a 2 ) =: a (ef' 2 v^ (27) 

one obtains after a short calculation the relation 

a = T^oe^ 2 (28) 



Reynolds' DNS data suggest Co ~ 7 |15|] , hence s 2 = 3 yields ao ~ 6 
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agreement with Bodenschatz's measurements |J. Moreover, for the Gaussian 
random variable Y = log ^ studied in eq. (|T6| ) one obtains the relation 



Y = --\og^---s 2 (29) 
2 6 (e) 3 V 7 

independent of the value of a and C . 

However, in general there are also other possible dynamical realizations 
of a Sawford model with fluctuating parameters than eq. (^]). For example, 
we may keep 7 constant in eq. ( p^) and let only a in eq. fluctuate. This 
leads to /3 ~ e -2 . Or, we may keep a constant and let only 7 fluctuate. This 
leads to f3 ~ e 1 / 2 . All these cases have in common that (3 ~ e K , where k is 
some power. Of course, for a log-normal distribution the power k does not 
change the functional form of the log-normal distribution, which once again 
is a hint for the physical relevance of a log-normally distributed (3. While all 
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the above possibilities lead to the same marginal distribution p(a) given by 
eq. they can be dynamically distinguished by looking at the acceleration 
autocorrelation function, which exhibits different behaviour for the various 
cases. 

To summarize, we have shown that the measured densities in Lagrangian 
turbulence experiments are very well described by log-normal superstatistics. 
This superstatistics is easily dynamically realized by considering a log-normal 
generalization of the model previously studied in 0. Only one parameter s 
is fitted to obtain excellent agreement with the experimentally measured 
distributions and the DNS data. Log-normal superstatistics differs from 
X 2 -superstatistics, i.e. ordinary Tsallis statistics, but for moderately large 
accelerations Tsallis statistics is often a good approximation. A possible hy- 
pothesis is that the variance parameter s 2 converges for R\ — > oo and is given 
by s 2 = 3. In that case no free parameters would be left, and a universal 
Lagrangian small scale statistics would arise. 
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Fig.l Acceleration distributions as measured by Bodenschatz et al. [|T], || 
for Reynolds number R\ = 690 and 970, and the log-normal superstatistics 
distribution p(a) given by eq. (§ with s 2 = 3.0. 
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